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ABSTRACT 

We explore nonlocally modified models of gravity, inspired by quantum loop 
corrections, as a mechanism for explaining current cosmic acceleration. These 
theories enjoy two major advantages: they allow a delayed response to cosmic 
events, here the transition from radiation to matter dominance, and they 
avoid the usual level of fine tuning; instead, emulating Dirac's dictum, the 
required large numbers come from the large time scales involved. Their solar 
system effects are safely negligible, and they may even prove useful to the 
black hole information problem. 
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1 Introduction 



A variety of complementary data sets pQ have led to general agreement that 
the universe is accelerating as if it had critical density, comprised of about 
30% matter and 70% cosmological constant [2j. There is, however, no current 
compelling explanation for either the smallness of A, or for its recent domi- 
nance in cosmological history [3]. Two existing classes of models, scalars jl] 
and u f(R)" modifications of gravity [5], can be arranged to reproduce the 
observed (or any other) expansion history [61 El E] . However, neither has 
an underlying rationale nor do they avoid fine tuning [9]. Quantum scalar 
effects, depending on a very small mass have also been proposed [TU] . 

In this work, we account for the current phase of acceleration through 
nonlocal additions to general relativity. Such corrections arise naturally as 
quantum loop effects and have of course been studied, though in other con- 
texts [6J [HJ [12] • As we will see, even the simple models we explore here can 
both generate large numbers without major fine tuning and deliver a delayed 
response to cosmic transitions, in particular to that from radiation to mat- 
ter dominance at z ~ 2300. We will neither attempt to derive our models 
from loop corrections nor to survey generic candidates here. Instead, we will 
show that natural nonlocal operators such as the inverse d'Alembertian can 
explain the time lag between z ~ 2300 and the onset of acceleration at red- 
shift z ~ 0.7, without recourse to large parameters. Large numbers come 
in our models precisely from the long time lags themselves, a mechanism 
reminiscent of some old ideas of Dirac. 



2 Nonlocal Triggers 



For simplicity, we deal with homogeneous, isotropic and spatially flat geome- 
tries 

ds 2 = —dt 2 + a 2 {t)dx ■ dx , (1) 
These correspond the following Hubble and deceleration parameters 

,,, s h , . oxl H , . 

H(.t) S - , ,( t ). -- = -!--, (2) 

and to Ricci scalai0 

R = 6(1 - q)H 2 . (3) 



iQur conventions are R = g^R^ and R^ = dpT^ - d„r% + TP T" - TP^T" 
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For much of cosmic history a(t) grows as a power of time 



a(t)~t s =► H(t) = - , q(t) = -. (4) 

t s 

Perfect radiation dominance corresponds to s — |, and perfect matter domi- 
nance to s = |. The Ricci scalar of course vanishes for s = | and is positive 
for s — |. It is the lowest dimension curvature invariant, and the only simple 
curvature invariant to vanish at finite s, so we concentrate here on i?-based 
models. 

We seek the inverse of some differential operator to provide the required 
time lag between the transition from radiation dominance to matter domi- 
nance at t eq ~ 10 5 years. The simplest choice is the scalar wave operator, 
suggested also by the fact that, for our background ([I]), dynamical gravitons 
obey the scalar wave equation [13] with 

Acting on any function of time fit), its retarded inverse reduces to simple 
integrations: 

= $[/](*) = - H dt '^) Jl ' dt"a\t")f{t") . (6) 

If we make the simplifying (and numerically justified) assumption that 
the power changes from s = \ to some other value at t = t eq , the integrals in 
([6]) are easily carried out for our choice of / = R 

g[R)(t) = _6^-i)( ln( M l_ _J_^-i| 

! - (3s -1) I H eq > 3s -1 3s-lW^ \ 



For the matter dominance value of s = 2/3, and at the present time of 
to ~ 10 10 years, this yields 

g[R](t ) _ a ~-14.0. (8) 

If we think of correcting the field equations by this term (apart from small 
additions that enforce conservation, and whose form we will shortly exhibit) 
times the Einstein tensor, this result already illustrates how nonlocality al- 
lows simple time evolution to generate large numbers without fine tuning. 
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Much larger values can be obtained through other operators, for example, 
the Paneitz operator arising in the context of conformal anomalies [H]. When 
specialized to our geometry ([I]) it takes the form 

-j= A F S tf + 2D^ - \,TR)D. ^ ^(444) ■ <9) 

One gets about 10 6 from the dimensionless combination of the inverse of this 
operator acting on R 2 . 



3 Specific Models 



Here we evaluate the consequences of the simplest alteration of the Einstein 
action, 

AC^^R^xf(Q[R]) . (10) 



16ttG 

One could modify the cosmological term in a similar way, but that turns out 
to require fine tuning to delay the onset of acceleration sufficiently 

Naively varying a nonlocal action such as ( flOl) would result in advanced 
Green's functions as well as the retarded ones (jEJ) we desire. However, be- 
cause conservation only depends on the Green's function being the inverse 
of a differential operator, one gets causal and conserved equations by simply 
replacing the advanced Green's functions by the retarded ones [12] J§ The 
resulting correction to the Einstein tensor is 



AG 



G>„ + g^D - D^D U ] ^f(Q[R]) +G[Rf'(G[R])] j 

d p (g[R})dJg[Rf(g[R}) 



As promised, it takes the form of a nonlocal distortion of the Einstein tensor, 
plus additional terms which enforce the Bianchi identity for any g^ u . The 
additional terms involve derivatives, so they are typically small when f(x) 
varies slowly Note also that, except for the very special case of f(x) = —x, no 



2 To derive causal and conserved field equations from quantum field theory one uses the 
Schwinger-Keldysh formalism [15] . This will generally result in dependence upon the real 
part of the propagator, as well as the retarded Green's function, which, if anything, may 
lead to even stronger effects than those we consider. 
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model of this form can be obtained from integrating out a scalar. Whatever 
these models' origin, then, they are not scalar-tensor gravities in disguise. 

Now note from (j7]) that Q[R]{t) is small for a long time after the on- 
set of matter dominance. During this period we may think of AG^ V as 
a perturbation of the stress tensor source, with AGoo = —SnGAp and 
g^AGij = —247rGAp. Our corrections will tend to induce acceleration if 
evolution during matter domination carries us to the point where 

AGqo + g ij AGij = -8ttG(A P + 3Ap) > QqH 2 = \-^- (12) 

Naturally, once our corrections exceed the Einstein range, they are no longer 
perturbations and numerical integration of the field equations is required. 
One illustrative class of models has 

f( x ) = Ce~i kx . (13) 

The resulting modification AG^ gives 

AGoo + g ij AGij ~ ^ • ± x c(l + h) (2 - 3k) (^-)\ (14) 

Note that the right-hand side is positive for k in the range — | < k < +|; 
actually, the range < k < +| is needed to make the correction term grow. 
Our results do depend on two dimensionless coupling constants, C and k, 
but neither need be very different from unity to provide a suitable delay for 
the onset of acceleration. For example, taking k = .1 and C = .2 would 
result in about the right onset time, in accord with the usual meaning of no 
fine tuning as involving parameters ~ 1. 

As stated earlier, it is possible to construct the scalar potential V(<p) to 
support an arbitrary expansion history a(t) obeying H > [BJ |7], and a 
similar construction exists for f(R) theories [SI [S]. The same possibility is 
of course present in our models, and indeed a procedure has recently been 
worked out for reconstructing the nonlocal distortion function f(x) which 
would support an arbitrary expansion history [16J. Hence there are certainly 
models of the type (ITUj) that fit the supernova data. Nor must one even resort 
to exotic choices of f(x). As might have been guessed from viewing these 
models as effective nonlocal distortions of Newton's constant, quiescence at 
recombination requires that f(x) be small for x near zero, whereas obtaining 
de Sitter expansion at asymptotically late times requires that f(x) approach 
— 1 from above for large, negative x. The onset of acceleration is controlled 
by the range of x at which f(x) becomes of order —1. 
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4 Conclusions 



We have explored the cosmological effects of some very simple nonlocally 
modified Einstein models inspired by loop corrections. Since their actual 
derivations from realistic quantum effects are likely to require nonpertur- 
bative summations, we regard them as purely phenomenological for now. 
Their two - equally important - main virtues are (unlike local variants): 
they naturally incorporate a delayed response to the transition from radia- 
tion to matter dominance, yet avoid major fine tuning. There are of course 
many other open questions raised by the present proposal, such as finding 
optimal candidate actions while ensuring that nonlocality has no negative 
unintended consequences. Some apparent worries, such as (unwanted) so- 
lar system effects, are easily allayed. There, G[R] ~ GM/(c 2 r) is a small 
number. Although a single power of Q[R] is observable — and constrains 
Brans-Dicke theory tightly [17] - higher powers, such as occur here, are 
negligible. 

It should also be mentioned that nonlocality may have a positive use in 
the black hole information problem [IB] : The infalling matter that creates 
or accretes to a black hole is imprinted on the external geometry through its 
stress tensor. Nonlocal dependence on the Einstein tensor will retain that 
information; while T^ u does not completely subsume the matter's internal 
structure, it is a significant repository thereof; furthermore, Q is singular on 
null surfaces such as the event horizon. 
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